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Abstract. Using four- wave mixing in a hot atomic vapor, we generate a pair of 
entangled twin beams in the microsecond pulsed regime near the Dl line of 85 Rb, 
making it compatible with commonly used quantum memory techniques. The beams 
are generated in the bright and vacuum-squeezed regimes, requiring two separate 
methods of analysis, without and with local oscillators, respectively. We report a noise 
reduction of up to 3.8 ±0.2 dB below the standard quantum limit in the pulsed regime 
and a level of entanglement that violates an Einstein-Podolsky-Rosen inequality. 
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1. Introduction 

Non-classical light is a key resource for quantum information and communication pQ 
as well as precision measurement [2]. Just as continuous- wave (CW) squeezing can be 
useful for noise reduction in gravitational wave detection [3j [4] , a robust source of pulsed 
squeezed light would be useful in the domain of quantum communication [lj. Indeed, the 
ability to store continuous variable entanglement in two separated quantum memories 
would allow for the entanglement of two macroscopic atomic ensembles [5j, similar to 
what has been achieved in the single photon regime [6j [7j. An experimental issue is 
the availability of an entangled light source compatible with a given quantum memory. 
The compatibility is often limited due to optical wavelength and bandwidth mismatch 
between them. 

In this paper we focus on generating quantum states of light suitable for Rb 
atomic memories such as those based on Electromagnetically Induced Transparency 
or Gradient Echo Memory [8J El HU] by using a Rb vapor near resonance as a non-linear 
medium. Four-wave-mixing (4WM) in Rb vapor has been shown to produce strongly 
correlated multi-spatial mode "twin beams" [HI [12l [T3l [HI [15]. Previous studies of 
4WM performed in the nanosecond pulsed regime have demonstrated squeezing [16j [IT] , 
but the inseparability of such states in the pulsed regime has not yet been established. 
Other techniques (e.g. Kerr effect in optical fibers) have successfully reached levels of 
entanglement that violate an Einstein- Podolsky- Rosen (EPR) inequality [18], and here 
we use an atomic medium to demonstrate up to 3.8 ± 0.2 dB of squeezing as well as 
EPR levels of entanglement between twin beams in the microsecond pulsed regime. In 
particular, we describe in detail two methods of analysis useful for characterizing bright 
beam intensity-difference squeezing and measuring entanglement between vacuum- 
squeezed twin beams, similar to [19], [20]. We do this by taking the Fourier transform of 
the noise time trace or integrating the photocurrent over the pulse window, respectively. 

2. Four- wave mixing and squeezed light 

Our experimental setup is based on the one described in [12]. We generate strongly 
correlated CW twin beams using four-wave mixing in a double-A configuration in a 
warm 85 Rb vapor. As illustrated in Fig. [TJ the four- wave mixing (4WM) process converts 
a pair of photons provided by a strong pump field into a pair of photons emitted into 
twin fields referred to as the probe and the conjugate. The pump beam is detuned 
approximately 800 MHz to the blue of the \5 2 S 1/2 ,F = 2) -> \5 2 P 1/2 ,F = 3) 
transition at 795 nm and focused down to a ^ diameter of 1.2 mm inside a 1.25 cm long 
85 Rb cell heated to 110°C. Taken individually, the probe and conjugate beams behave 
like thermal states exhibiting excess quadrature noise relative to the shot noise limit 
(SNL). When compared these twin fields are known to exhibit strong correlations [21J. 

To obtain bright twin beams, we stimulate the 4WM process by injecting a coherent 
state "seed" into the probe mode, which produces bright probe and conjugate beams 
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suitable for direct detection with photodiodes (Fig. [Tp). These bright beams are 
intensity-difference squeezed, and the squeezing is measured by directly detecting the 
probe and conjugate intensities with matched photodiodes with quantum efficiencies of 
95%. The photocurrents are subtracted and the intensity-difference signal is amplified 
and recorded. 

Although easy to implement, direct detection does not provide access to 
measurements of the phase quadrature, a necessary condition to show entanglement [22J. 
To measure these anti-correlations, a common approach is to perform a simultaneous 
homodyne detection of each beam (Fig. [T^) . The generalized quadrature Xq. can be 
defined as Xq. = -^(a^e l$i + ae~ %di \ where % can be p or c for the probe and conjugate 
beams respectively. The phase of each local oscillator (LO) can select the quadrature 
to be detected in each beam and accordingly, if the relative phase between the LOs is 
stable over the detection time, it is sufficient to sweep only one LO phase in order to 
fully characterize the state. 

When the signals obtained by the homodyne detections are subtracted, the noise 
power of the joint quadrature Xq = X$ p — X$ c where 6 = 9 P + 9 C is measured. In 
this case it is possible to observe squeezing provided that 9 = 9~ = 0. Similarly, 
we can measure the sum of the two homodyne signals in order to obtain the joint 
quadrature X$~ = X$ p + Xfj c which is squeezed if 9 = 9 + = tt/2. The observation 
of squeezing on both joint quadratures is sufficient to demonstrate entanglement [22J. 
The degree of entanglement can be quantified by the inseparability /, defined as 
the sum of the noise variances measured for X^ + and X^L normalized to the SNL : 
I = ((A!") 2 ) + ((AX+ ) 2 ). States are said to be entangled if J < 2 [231 EH . 




Figure 1. Four- wave mixing in 85 Rb vapor, (a) Energy levels, (b) Vacuum-squeezed 
twin beam generation. Probe and conjugate fields are coupled over a range of angles 
Acj) determined by the phase matching conditions, and selected by the local oscillators 
used for homodyne detection, (c) Bright twin beam generation with a coherent state 
seeding the probe field at </>. 
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3. Intensity-difference squeezing 

We apply these techniques to the pulsed regime for bright beams and vacuum-squeezed 
twin beams. We first measure intensity-difference squeezing between pulsed bright 
beams as sketched in Fig. [2| A coherent state at the probe frequency is obtained by 
double-passing an acousto-optic modulator (AOM) coupled to a 1.5 GHz radio frequency 
(RF) signal. The RF drive is then pulsed in time, producing a train of pulses for the 
input seed. While the applied RF pulses were square pulses, the optical pulses were 
bandwidth-limited by the 4WM process. High-frequency components of the pulses are 
not amplified by the process and have a negligible contribution to the noise measurement. 
These pulses are shot noise limited which is an important condition in order to observe 
intensity-difference squeezing. The pulse width is 2 /is and pulses are generated with a 
period of 10 /is. These bright pulses are directly detected with matched photodiodes and 
the resulting photocurrents are subtracted. The noise of the subtracted photocurrents 
therefore alternates between the intensity-difference noise of the squeezed light and the 
electronic noise floor of the detection system when the pulses are not present (Fig. [2]). 

The arrival of the probe lags the conjugate by approximately 10 ns due to the 
dispersion in the atomic vapor [25J. This difference in the arrival times generates 
a voltage spike from the RF amplifier, and the capacitance of the detecting circuit 
produces a large ringing effect that dominates the intensity difference signal at 
frequencies below 0.5 MHz. A delay line does not resolve this issue because the exact 
time delay between the pulses fluctuates and it is not possible to cancel the voltage spike 
to the required level. Although the noise power introduced by the ringing primarily 
dominates frequencies below 2 MHz, it is large enough to saturate the amplifying 
circuitry of the detection system. To suppress this ringing, we included a high-pass 
filter with a cutoff frequency of 300 kHz in the detection circuit. 

We detected the shot noise by pulsing a shot-noise-limited beam picked off from 
the master laser with the same duty cycle as the seed, splitting the pulse with a 50- 
50 beamsplitter, and directly detecting with the same detection technique and the 
same detection circuitry (see Fig. [2]). Given the detected intensity-difference noise, the 
intensity-difference squeezing was computed as indicated in Fig. [2J3. Each portion of the 
amplified photocurrent corresponding to the noise of the twin beam intensity difference 
(blue regions in Fig. J^) consisted of 200 values sampled over 2 /is, adequately sampling 
frequencies over the w 20 MHz bandwidth of the 4WM process. We computed the 
Fourier transform for each blue region in Fig. [2)3 and then averaged 10 3 power spectra 
to obtain the final power spectrum. 

The average pulsed intensity-difference squeezing spectrum is plotted alongside CW 
measurements in Fig. [3| The CW measurements for the shot noise and squeezing were 
taken without the high-pass filters pictured in Fig. [2] since the ringing effect discussed 
above does not apply to the CW case. As shown in Fig. [3j the intensity-difference 
noise power in the pulsed and CW regimes is similar above 3 MHz. At lower analysis 
frequencies, however, we observe a discrepancy between these two regimes in the form 




Figure 2. Experimental procedure to detect pulsed bright beam intensity- difference 
squeezing. PBS, polarizing beamsplitter; BS, 50-50 beamsplitter; HPF, high-pass filter, 
(a) Detection of pulsed bright beam intensity-difference noise, (b) Pulsed shot noise 
detection, (c) Analysis procedure to measure intensity-difference squeezing between 
bright beams. The pulse train of bright beams yields a subtracted photocurrent signal 
(red trace) that alternates between the squeezed intensity-difference signal (regions 
shaded blue) and the electronic noise floor (regions shaded green). The Fourier 
transform of the time trace in each blue region is computed and then averaged to 
obtain the final power spectrum of the intensity-difference noise. Traces in this figure 
are not experimental results and are included only to clarify the analysis procedure. 



of a large excess noise for the pulsed measurements. We suggest that the excess noise 
on the squeezing partially comes from the fact that the residual ringing at the detection 
stage described earlier is not completely filtered by the high-pass filter. At frequencies 
larger than 10 MHz, CW and pulsed squeezing vanishes and the noise level rises above 
SNL. This effect is due to a different propagation velocity of the probe and the conjugate 
inside the atomic medium and can be easily mitigated using an electronic delay line[26j. 

4. Vacuum-squeezed twin beams 

In order to demonstrate that the beams are entangled, we would like to generate pulsed 
vacuum-squeezed twin beams suitable for homodyne detection. To obtain these pulsed 
beams, we first generate CW vacuum-squeezed twin fields (Fig. [TJd) and direct the 
probe field through two AOMs that act as a switch. These AOMs are driven with 
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Figure 3. Intensity-difference squeezing obtained for the the continuous- wave (red) 
and pulsed (blue) regimes. The shot noise (black) and electronic noise (green) are 
plotted for reference. The data presented here are raw data without any corrections. 
When these data are corrected for the electronic noise, the value of the shot noise is 
constant from 100 kHz to 15 MHz. 

synchronously pulsed RF power to deflect the mode when the squeezed vacuum is to be 
turned off, and oriented such that the acoustic waves propagate at right angles to each 
other to increase the extinction of the mode(Fig.[3^). The conjugate beam remains CW. 
The probe and conjugate are then detected with separate homodyne detectors with CW 
LOs and the quadrature values recorded. It should be noted that pulsing the 4WM 
pump beam is not an option to produce pulsed vacuum-squeezed twin beams since the 
4WM process typically requires several [is of pumping time to reach steady state, which 
is on the time scale of interest for the pulses. 

Nonzero diffracted orders emerging from the AOMs are blocked with irises placed 
immediately after each AOM. We choose to detect the non-diffracted order since it suffers 
less loss and less mode distortion than the deflected beams. Preserving the spatial mode 
quality of the probe is critical since any mode mismatch between the probe field and the 
probe LO used in the homodyne detection contributes to an additional effective loss. The 
disadvantage of this technique is that extinction of the probe is limited by the diffraction 
efficiency that can be reached by the AOMs. Orienting the AOMs perpendicularly to 
one another allowed us to achieve 97% suppression of the probe mode, when activated, 
and a transmission greater than 95%, when not activated. 

To measure the noise of the phase-dependent joint quadratures X$ and Xg~, we 
pulse the probe field 10 4 times over 100 ms with a period of 10 /xs (Fig. ^p) as we 
sweep the phase of the probe LO by n. This technique assumes a certain degree of 
phase stability of the probe LO over the measurement time p~9j EOj [27] . To verify 
that the phase excursion between the probe and conjugate is negligible over the relevant 




Figure 4. (a) Detection of vacuum- squeezed twin beams, (b) Procedure to measure 
joint quadratures, (i) The CW homodyne detection of the conjugate (blue trace) 
samples the conjugate generalized quadrature Xq c . The probe homodyne detection 
(green trace) alternates between detecting the probe generalized quadrature X v Qv (blue 
regions) and vacuum (green regions). Probe and conjugate traces are recorded. An 
offset time At is introduced when the data are analyzed in order to maximize the 
squeezing, (ii) The difference signal (red trace) between the probe and conjugate 
homodyne detection alternates between the joint quadrature Xq (blue regions) and 
the sum of the conjugate excess noise and the vacuum noise (green regions). Similar 
operations are needed to measure X#~. (iii) We spectrally filter Xq by multiplying the 
blue region with the window function of Eq. [I] (iv) The values of Xq are obtained 
after integration over the blue regions. At the end of the measurement, the shutter 
controlling the probe and conjugate fields is closed and the shot noise (red regions) 
is measured analogously. Traces in these figures are not experimental results and are 
included only to clarify the analysis procedure. 
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timescale, we seed the twin beams with a coherent state and detect the resulting classical 
interference with the LOs. Over periods of 100 ms, we typically obtain a phase standard 
deviation on the order of 1 degree. 

When the probe is turned on, the homodyne difference signal yields a detection of 
Xq and the sum of Xq~. When the probe is off, the probe homodyne detector measures 
shot noise and the conjugate homodyne detector measures the uncorrelated generalized 
quadrature Xq c . The difference and sum of these two contributions yields the same noise 
power in excess of the two-beam SNL. At the end of the 100 ms sequence, the twin fields 
are blocked with a mechanical shutter and we record 10 ms of two-beam vacuum shot 
noise. 

As discussed elsewhere [HJ [21], the 4WM process exhibits a frequency-dependent 
squeezing spectrum over a bandwidth of approximately 20 MHz. In order to study the 
phase-dependence of the noise power of the joint quadratures, we wish to confine our 
attention to a limited bandwidth where the squeezing is maximized. To achieve this 
goal, we multiply each portion of the time trace corresponding to the measurement of 
the joint quadratures (blue regions in Fig. |i]>iii) by the window function : 

W(t) = — ^= cos(uj (t - t ))e^^e(r/2 - \t - t |) (1) 

0(t) refers to the Heaviside step function, r the probe pulse width, and i the midpoint 
of the probe pulse in time. According to the convolution theorem, multiplying the noise 
time trace by W(t) allows us to study a band of frequencies in the squeezing spectrum 
centered around ujq/2tt with a ^ spectral width of a -1 . CW measurements of the 
squeezing reveal that the squeezing is maximized over a bandwidth of approximately 300 
kHz centered around 750 kHz. Accordingly, we selected a value of 0.5 [is for the window 
width a. We plot 10 4 recorded values of Xq versus 9 in Fig. [5^i (red points) alongside 
a sampling of the twin beam shot noise (black points). The same data were recorded 
for Xq~ but are not shown for the sake of clarity. To quantify the level of squeezing, 
the red points in Fig. [5^l are binned over 100 equally-spaced phases of the generalized 
quadrature 9. The variance of the Xq and Xq~ values in each bin was computed and 
plotted against the average joint quadrature phase 9 of each bin in Fig. ^jp. The SNL 
is obtained by computing the variance of the black points in Fig. [5^. From the figure, 
we observe 3.8 ± 0.2 dB of squeezing for both Xq_ and Xq~+. The uncertainty of 0.2 dB 
is obtained by calculating the standard deviation of the noise power over a phase range 
of 7r/20 around the phase 9 = n/2. The measured inseparability is I = 0.83 ± 0.04 < 2, 
which is sufficient to conclude that the beams are entangled. 

A more stringent criterion is commonly used to describe a higher degree of entan- 
glement, known as EPR entanglement. This criterion quantifies the knowledge about 
the state of the probe beam by measuring the conjugate. Although conditional variances 
are needed to quantify EPR entanglement, for Gaussian twin beams it has been shown 
that the product of the conditional variances is bounded by 4 ^(AX^_) 2 ^ ^(AXj^) 2 ^ and 
therefore a value smaller than 1 for this product establishes EPR entanglement. From 
the data shown on Fig. [5j3, we obtain an upper bound for the product of conditional 
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Figure 5. Detected vacuum squeezing, (a) Normalized values for the joint quadrature 
Xq of the twin beams (red) as a function of the joint quadrature phase and two- 
beam vacuum shot noise calibration (black), (b) The measured quadrature values are 
binned over 100 equally spaced phases of 0, and the variance (AX^) 2 of each bin is 
plotted against the 100 corresponding average values of 0. The black line indicates the 
shot noise limit (SNL). 

variances of 0.69 ± 0.07 < 1, which is sufficient to verify the presence of EPR levels of 
entanglement. 

5. Conclusion 

In this paper we demonstrate a source of pulsed bipartite entangled beams based on 
four- wave-mixing in a hot atomic vapor. Substantial improvement in the degree of 
correlation (3.8 ± 0.2 dB below the SNL) and bandwidth (MHz) compared to previous 
work with an atomic vapor is reported. Inseparability and EPR levels of entanglement 
are demonstrated with this technique in the pulsed regime. While comparable results 
have been reported using the Kerr effect in an optical fiber at 1.5/iin [18], in this paper 
we demonstrate pulsed beams that match the bandwidth and the wavelength of atomic 
quantum memory techniques involving Rb vapor. Finally, as 4WM is known to exhibit 
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multi-spatial-mode properties [12], we think that this source will be of great interest 
for spatially multimode atomic memories [28] . 
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